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Abstract 

A global extension theorem is established for isotropic singularities in polytropic 
perfect fluid Bianchi space-times. When an extension is possible, the limiting be- 
haviour of the physical space-time near the singularity is analysed. 
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1 Introduction 



In an earlier article [T] we proved a local conformal extension theorem, estab- 
lishing curvature conditions along an incomplete conformal geodesic sufficient 
to permit an extension of the conformal structure. The proof was phrased 
in the tractor formalism. As was noted at the end of [1], one would like to 
generalise this result to a global extension theorem, extending simultaneously 
through all singular points. In this article we analyse perfect fluid Bianchi 
space-times for isotropic singularities. We find that the fiuid fiow-lines are 
conformal geodesies, which allows us to apply the local extension theorem of 
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[T] to any one of them. Then the spatial homogeneity makes it possible to 
obtain a global extension result. 

Isotropic singularities in perfect fluid space-times were previously studied by 
using the congruence of flow lines. In particular, [2] analysed the limiting 
behaviour of the physical space-time near the singularity when the latter is 
isotropic. It will be shown in Section[5]how these limits can be directly obtained 
from boundedness conditions on the tractor curvature. 

We use the 1 + 3 decomposition with respect to a time-like congruence, which 
we outline in the next section. Our setup will follow that of Pl5|6] , so that for 

this paper (unlike [1]) we will adopt the signature convention ( h ++), as 

used in this literature. 

We make use of conformal densities and make some use of the tractor for- 
malism, though most formulae are tensorial. For an introduction to tractors 
and the notation the reader is refered to [1]. In abstract indices, the confor- 
mal, physical and unphysical metric are written gab, dab, dab respectively, so 
that gab = ^^Qab- To simplify notation we write all equations in the untilded, 
unphysical variables (V,(7). However when particular physical assumptions 
simplify the equations, we will highlight this by using the physical variables 



2 l-|-3 decomposition 

We review those parts of the formalism of 1 + 3-decomposition with respect to 
a time-like congruence that we shall need, indicating the conformal properties 
of the familiar kinematic quantities associated with the congruence. Our main 
interest is in the spatially-homogeneous perfect-fluid cosmologies, when we use 
the fluid congruence, which is geodesic and twist-free in the physical metric. 
Since we shall only rescale with functions of t, it has these properties in the 
unphysical metric too, and it also turns out to be a congruence of conformal 
geodesies. 

Given the time-like congruence with velocity field v°' we define the velocity 
scale V := \gabV°'v^\^ ^ in the notation of [1]; recall this means that 

u is a conformal weight one scalar, a section of the line bundle r(e[l]) of 
these. Then there is an associated unit velocity u"" := ^ G r(e"[— 1]) and 
acceleration := u"-Vau'' G r(£:^[— 2]), respectively conformal weight -1 and 
-2 vector fields. Given the application we have in mind, we will set the twist 
of the congruence to zero from the outset. 

We have projection operators: = —UaU^ onto the flow lines and hj' = 
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5 a —Uj" orthogonal to them; thus U"^ = U, = h and Uoh = hoU = 0. 
The operator /i"^ := h^"'Ji \ — \h°-^hcd projects out the spatial, symmetric, 
trace-free part of a tensor. All three operators have conformal weight zero and 
are used to decompose tensors into components along and orthogonal to the 
flow. 

With the twist-free assumption, the decomposition of the covariant derivative 
VaMfc takes the simplified form 

VaUfe = -UaAb + "^dhab + (Tab € T{eab[^]), (1) 

in terms of the kinematic quantities: the acceleration G F (£:''[— 2]), expan- 
sion 6 G r(£:[— 1]), and shear aab G r(£:(ab)[l]). 

We define the magnitude = \(Jab<^°'^- The flow is hypersurface-orthogonal, 
with induced metric hab and second fundamental form Xab = \dhab + CTab on 
these hypersurfaces. 

We recall (pQ) that a conformal density is homogeneous if the change under a 
gauge transformation to a general Weyl connection V = V -|- has no additive 
correction in the one-form he which determines the change in connection. We 
can express the inhomogeneous conformal density 9 in terms of an average 
length scale L G r(£:[l]) by 

e = ^Y^^e = e + -i{h,u)eT{e[-i\), (2) 

Ij 

and then L is a homogeneous density. 

The Einstein tensor Gab is related to the Schouten tensor Pab by 

Gab '■= '^{Pab — PcdO^'^gab) 

and in the physical space-time is equated to (a multiple of) the energy mo- 
mentum tensor Tab- Under connection translation V = V + 6 the Schouten 
tensor transforms as 

Pab - Pab = ^abb " bah + ^gabg'Xbd- (3) 

In this article we will only use Levi-Civita connections, so that the 1-forms ba 
are always closed and the Schouten and Einstein tensors are always symmetric. 
In general, therefore, we deflne the geometric matter variables by decomposing 
the Einstein tensor Gab according to 

Gab = fJ^UaUb + QaUb + UaQb + phab + T^ab, (4) 
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where /i and p are scalars, qa is a vector orthogonal to and iiah is a trace- 
free, symmetric tensor orthogonal to u"". (In [2] the geometric matter variables 
were denoted by A, B, S^, T,ab-) In the physical space-time, the corresponding 
matter variables are tilded: fi is the energy density relative to tt", p is the 
isotropic pressure, and and Tfab are the momentum density and the trace- 
free anisotropic pressure, both zero for a perfect fluid. The trace-free Ricci 
tensor ^^^Sab of the hypersurfaces can be written as 

11 2 

^^^Sab = Eab + -TXab - -^CTab + (^a^^cb " -'^'^Kb- (5) 



where Eab is the electric part of the Weyl tensor. 

Finally in this section, for a perfect fluid the conservation equation implies 
-3^- (6) 



which can be solved for a polytropic perfect fluid, when p = (7 — l)p, to give 
p. = fi,L-^\ (7) 

where p.^ is constant along the flow. 



3 Bianchi space-times 



A space-time is said to be spatially homogeneous if it admits a 3-dimensional 
vector space of Killing vector fields (KVFs) transitive on space-like 3-surfaces. 
The KVFs form a Lie algebra with basis ^a' say, with a' = 1,2,3, and the 
associated Lie group G3 is the corresponding group of isometrics. The Bianchi 
models are classified in terms of the structure constants C^^,^, of the Lie 

algebra, where [^a',^/?'] = ^,^,^^1 (see e.g. |6H7] )and one then finds nine 
canonical forms, the Bianchi types. 

We shall restrict to non-tilted perfect fluid cosmologies, for which the velocity 
u of the fluid (a unit vector for the physical metric) is orthogonal to these 
3-surfaces, which can be labelled by proper-time t along the fluid flow. Co- 
moving coordinates x- , with a' = 1, 2, 3, (coordinate indices underlined) can 
be introduced in a standard way for each Bianchi type. We denote the Jacobi 
fields obtained from these standard comoving coordinates by 77^/ = d/dx- . 
From spatial homogeneity, we have C^u = 0, that is the velocity is invariant 
under the action of the KVFs. To construct a group-invariant orthonormal 
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frame at one time, choose an orthonormal frame {e^}, a = 0,1,2,3, at one 
point on an orbit with Cq = "U* and Lie drag it around the 3-surface, so that 
C^e\ = 0. Then C^g{ea, (ip) = so that the frame is orthonormal everywhere 
on the orbit. 

For the propagation along the congruence, note that £[ji,|]e^ = so that 
necessarily 

= M/{t)i^p, (8) 

for some M^^{t), which at this stage is freely specifiable. 

In the next section we shall see that the fluid flow lines are conformal geodesies, 
when it is natural to choose a group-invariant frame which is Weyl-propagated, 
that is 

V.el = 0, 

where V is the Weyl connection defined from the conformal geodesic equation 
and is tangent to the conformally-parametrised conformal geodesies. We 
shall take u = Cq. Then 

where 9ij is the unphysical second fundamental form which, by the spatial 
homogeneity, depends only on time. The connection coefficients T^^^^ of the 
connection V in the chosen frame are defined by 

Ve,e^ = r%e;. (9) 



4 The fluid flow-lines as conformal geodesies 

In the physical space-time the conformal geodesic equations can be written in 
a third-order form (given in [1] , but remembering that the signature here is op- 
posite to the one used there): if the unit velocity vector is with acceleration 
y4* then 

DA' + {u''bAk)u' = -P/u^ - {PjkU^u'')u\ (10) 
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where D = u^Vi. The conformally-parametrised tangent vector is = g^w* 
where q is found from 

2^ = -P^^u^u^ + \a^A\ (11) 



and the one- form 6j, associated with the congruence and used in Weyl propa- 
gation, is obtained as 

bn 

hi = -Ai + (12) 



We set A = -2^ = h^i . 

q 1 

In a (non-tilted) perfect fluid Bianchi space-time the matter flow lines are 
space-time geodesies and the Schouten tensor has the form 

Pzj = (^/i + UiUj + ^fihij. (13) 



It is then straightforward to check that (fT0|) - (fT2|) can be solved with Ai = 
and bi = —Xui. Thus the flow lines of a perfect fluid in a Bianchi space-time 
can be reparametrised as conformal geodesies. The 1-form bi is exact and hence 
V = V + 6, the Weyl connection associated with the congruence, is actually 
the Levi-Civita connection preserving the metric Qab = ^'"^Qab for Q = q'^. In 
particular therefore the Weyl-propagated tetrad is parallelly-propagated in the 
unphysical metric. The conformal parameter r along the conformal geodesies 
is given by 




up to additive constant, where t is proper-time in the physical metric. In 
the rescaled space-time, the conformally-parametrised tangent vector = 
g^M* is a unit vector, so that it is in fact affinely parametrised, and r is the 
corresponding proper time in the unphysical metric. 

For a perfect fluid with a polytropic equation of state and polytropic index 7, 
(fTTl) takes the form 

.. 1/1. lA 37-1. 

^ = -2l3^+2^j'^ = "^2-^^' ^^^^ 
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with the overdot for d/dt. The freedom in the choice of solution to cor- 
responds to the freedom in fractional- hnear transformation of the conformal 
parameter. 

Suppose that, into the past along each worldline, we reach the singularity as 
t — ^> 0. The volume goes to zero in the physical metric, so for the volume to 
be nonzero in the rescaled space-time we require g — >^ in that limit. 

We next assume that there is a choice of q such that r does not have infinitely 
many poles before the singularity is reached, and is not automatically infinite 
at the singularity. As discussed in [1] , this is an assumption on the singularity 
in the physical space-time which is logically prior to that of bounded tractor 
curvature. Thus from now on we assume that q vanishes at t = but not in 
some initial interval (0,ti), and that the integral for r 

r = j -dt, (15) 





is finite for t < ti. 

From the definition of A, we find the expression q = goe"^/'^'^*, so that A — > 
— oo as t — >■ 0. Furthermore \q^ = —2qDq diverges near the singularity. To 
see this, suppose conversely that, on an initial interval (0,e), Dq"^ = —\q^ is 
bounded, by C say, then 

q^{t) < Ct so that j q-^{t)dt > J ^ on (0,e), (16) 



but this would violate our assumption on r after (IT^ . These limits are needed 
in Section O 



5 The Global Extension Theorem and limits at the singularity 

In [1] the following local extension theorem, for singularities to the future, was 
proven: 

Theorem 5.1 Let 7 : [0,Tp) ^ M be the final segment of an incomplete 
conformal geodesic in (M, g), such that b is bounded in [0, Tp). Let W C M be 
a neighbourhood o/7[0,ri?) in which the strong causality condition holds. Let 
{e^} be a Weyl propagated orthonormal frame along 7 with associated tractor 
frame Ss- 
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i) Suppose the tractor curvature TZ and its derivative V^TZ along'') have hounded 
norms on 7 with respect to {ep} and Eq. Then there exists a neighbourhood U 
o/7[0, Tp) with U <ZW and a diffeomorphism V (Z H'^ ^ U. 

a) Suppose that in U the tractor curvature norms H?^*^'^"^"'^-' || and H?^*^^''"'"'"^'' || 
are bounded. Then there exists a general Weyl connection V associated to [g] 
and a conformally related metric gij such that there exists U* D U with a 
C^-extension of {U,gij\u) into {U*,gij\u*)- 

Hi) The Riemann curvature of gij is C^~^ . 

Thus the conformal structure (M, g) is locally extendible. 

Definitions of tractor curvature and tlie norms used liere are given in [Tj, 
but tlie terms can be understood as follows: the components of the tractor 
curvature TZ are the Weyl tensor and the Cotton- York tensor; boundedness in 
the tractor frame is equivalent to boundedness in the Weyl-propagated tetrad; 
the norms are Riemannian sums-of-squares in the specified basis; and ||7^*-'''^^'' || 
and ||7^(^'^+l)|| are norms of {k + l)-fold derivatives of TZ and their derivatives 
(once) along u. 

In this article, we are considering initial rather than final singularities and 
hence we use initial segments (0, ri]. In [1] it was remarked that the extension 
theorem was only local because one could not guarantee that the set O, the 
neighbourhood of 7 used for the extension process, included more than one 
singular point. For Bianchi space-times, the spatial homogeneity ensures that 
all the conformal geodesies in the perfect fluid congruence are equivalent. 
Hence conditions satisfied by the chosen central curve 7 are satisfied by all its 
neighbours. In particular, this means the work underlying {€) of Theorem 15.11 
is avoided since the usual coordinatisation of Bianchi space-times provides a 
good coordinate system. The local extension theorem above will therefore give 
a global extension theorem when applied to perfect fluid Bianchi space-times. 

The version of the Whitney Extension Theorem used in pQ was for sets in 
which satisfied Whitney's Property P. We needed bounds on the derivatives 
of order k + 1 oi the function /, which it is desired to extend, in order to obtain 
appropriately uniform bounds on derivatives of / of order /c; then the theorem 
provides a C'^-extension F of /. In the spatially-homogeneous context, we only 
ever consider functions of the one variable t and so we have a simpler, in fact 
elementary, extension theorem: 

If f & C'^((0, Ti]) with all derivatives bounded and having limits at zero, then f 
admits an extension F to 00, ri]) which is in fact analytic on (— oo,0). 
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This is obtained from the polynomial 

in t < 0. Thus we may assume a lower degree of differentiability here. 

Finally the assumption of strong causality was made in [1] in order to prove 
Property V, so again we don't need to make it in this context. 

A global extension theorem with minimal conditions on the tractor curvature 
can be stated as follows: 

Theorem 5.2 Suppose that we are given a perfect fluid Bianchi cosmology 
with an initial singularity, and a choice of q and an initial interval for which 
the conformal parameter t is finite at the initial singularity. Suppose further 
that along one, and hence any, fluid flow-line the tractor curvature TZ and 
its derivative VuTZ along u are hounded in a b-propagated frame. Then there 
exists a conformal rescaling such that the rescaled metric gij and its inverse 
are hounded at the initial singularity and hence extendihle. Furthermore the 
unphysical quantities 

^'^Rab (17) 

defining the curvature and connection, are all hounded, with limits at the sin- 
gularity. 

Proof: The 1-form b is exact and hence the general Weyl gauge and the unit 
velocity gauge (unphysical gauge) coincide: V + 6 = V = V. The 6-propagated 
frame is parallelly propagated with respect to V, the flow lines are metric 
geodesies of gab and the vectors ria' are Jacobi fields. By Theorem 4.1 and 
Lemma 4.7 of [1], boundedness of TZ implies boundedness of the metric in the 
coordinate basis and the kinematic quantities in the Weyl-propagated basis, 
and that these quantities have finite limits at the singularity. Theorem 3.4 of 
[1] now gives boundedness of the inverse metric in the coordinate basis, and 
consequently of the matrix relating the Weyl-propagated basis to the Jacobi 
fields, and its inverse. These quantities also have limits at the singularity. 

The Weyl tensor is part of the tractor curvature and so is bounded by as- 
sumption. For the Schouten tensor (equivalently, the Ricci tensor) we recall 
the definition 

Yabc = 2V[a-P6]c, 

of the Cotton- York tensor, which is bounded by the assumption on TZ. With re- 
sults in hand, it is straightforward to integrate this to find that Pab is bounded. 
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with a limit at the singularity. Thus all of the unphysical Riemann tensor is 
bounded at the singularity. Now Lemma 4.11 of [1] with j = 1 implies bound- 
edness of V^e^ for all choices of rj, and therefore also of Ve^e^ (in fact one 
obtains boundedness of the tractor derivatives of the tractor frame). This gives 
boundedness of the connection coefficients r[^^ with a limit at the singularity. 

Finally, if VuTZ is also bounded then the Weyl curvature, which at this stage 
we know to be bounded, also has a limit at the singularity. From this follows 
that the 3-curvature of the hypersurfaces, ^^^Rab, is also bounded and has a 
limit at the initial singularity. □ 

The unphysical space-time is well behaved and bounded as long as the tractor 
curvature is. We use this fact and combine it with the rescaling transformations 
for curvature and the geometric matter variables to get some information on 
the unphysical space-time. 

Proposition 1 Suppose the conditions of Theorem 15.^ hold. Let ip = \q^jl 
and Ip = then the following limits hold as t 0: 

'-^l (18) 

2 u 1 5-^ , ^ 

e^+oo, --^77, ^-^0, 19 

/-^ r^abcd 

^^^^ > (20) 

- Sl^(T^f3 ^ 0, H^fs ^ 0. (21) 



Remark: The limits (fT9l - [20]) coincide with those derived for isotropic singu- 
larites in Theorem 3.2 and 3.3 of [2]. 

Proof: We use the boundedness of the curvature components to prove the 
limits above. The conformal geodesic equation (fTOj) gives us Pabv"" = and 
hence by (IT^ + |p = 0. Combining ba = —Xua = —^JUa with ([3D we have 



Pab = Pab - '^a'ipUb - IpV aUb - Ip'^UaUb - ^'ip'^gab, 
'^aUb = VaUb + baUb + Uah - QabbcU" = Bab - 24jUaUb - tpQab- 

Contracting the first equation with h"-^ = q'^h"-'^, respectively /i^^, and substi- 
tuting (fT3l) with the perfect fiuid condition on the left hand side, we obtain 



^ = P/-^e + ^V'', (22) 
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(23) 



where (!22l) uses (fTTI) . For the conformal time parameter to be finite at the 
boundary we have shown that ip is required to diverge. Since Pat is bounded, 
it follows that iVab is bounded and, from (l23l) that the unphysical shear vanishes 
as t — > 0. 

From fl22]) . we deduce that must diverge and then obtain the hmit f[T5]) . 
Next we note the transformation of the following kinematic quantities under 
rescaling: 



where now Q = q"^ and u°'Ta = —X . Using the boundedness of 6 and (!24l) 
and substituting for A we get the first two of (fT9!) as t — >• 0. For the third we 
use ([25D in the form cr^ = q'^a^ and the observation, made after (HM . that Xq"^ 
diverges. 

For the quotient of curvatures ([20]), we observe that CatcdC"-^^'^ = q'^CabcdC"^^'^ 
and CahcdC''^'"^ is bounded, while q^RabR"^ = g^(/i^ + 3p^) > 4(^^ and so 
diverges. Hence this ratio tends to zero. 

For a polytropic perfect fiuid in the physical space-time, ([7]) gives /i = fi^L~^'^ = 
fl^q~^'^ L'^'^ . Since 6 is bounded, L tends to a finite non-zero limit and there- 
fore so does q^'^fl. The range 7 < § and equation fl23|) would now imply that 
9?, are bounded. As discussed earlier, this would prevent us from reaching 
the singularity in finite r, so that this range in 7 has been ruled out by our as- 
sumptions (of course on physical grounds we are only interested in 1 < 7 < 2). 

If 7 = I then fiq^ has a finite non-zero limit. Then by ffTS!) Xq^ is finite too. 
Rewriting (1251) with the physical shear, Xq'^dab = h'^Pcd , we get a finite limit 
for dab- Similarly if 7 > | the physical shear must vanish at the singularity. 

To analyse the behaviour of the Weyl curvature we look at the tractor curva- 
ture. We assumed that V^TZ, TZ are bounded in the 6-propagated conformally 
orthonormal frame to be able to deduce the existence of a space-time extension 
through the initial singularity. It followed that Pij,Cijki and the Cotton- York 
tensor Yij^ are bounded. We have also shown that A and ip = Xq^ need to be 
unbounded. 




(24) 



(25) 
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In the physical space-time with perfect fluid we use ([6]) to deduce that the 
physical Cotton- York tensor is given by 

Yijk = 2W[iPj]k = -(/i + p)u[iaj]k. (26) 

From the transformation rule of the Cotton- York tensor we obtain that Yiji + 
^ij^^k = Yiji is bounded. In the Weyl-propagated frame we flnd 

= i)Ea^ - + p)(J„^, (27) 

l^pg%ie>'e\ = q'XH^, = ^H^,. (28) 

Equation ( !27l) relates the unphysical electric Weyl and shear tensor to the 
physical density and pressure. The boundedness of the left-hand-side of (1271) 
with (fTSj) gives the flrst of (12T|) ; the boundedness of the left-hand-side of (!28l) 
and the fact that ip diverges imply the second. 

From ([5]) at t = and using these limits we obtain 

^ 'Dab — i^ab + -T^ab — ^ab, 

2 37 

so that the trace-free 3-curvature of the hypersurfaces of homogeneity has a 
limit at the initial singularity, where it is determined by the initial electric 
Weyl tensor and vice versa. □ 

These limits for isotropic singularities were found in [2] by assuming the ex- 
istence of an initial isotropic singularity. We have derived them directly from 
the assumption of bounded tractor curvature and its derivative, and flniteness 
of the conformal time parameter, in the particular case of spatial homogeneity. 



6 Conclusion 

We have shown that perfect fluid flow lines in Bianchi space-times are con- 
formal geodesies, which has allowed us to apply the Extension Theorem 15.11 
Then spatial homogeneity has resulted in a Global Extension Theorem 15. 2[ 
where we can extend across the whole singularity at once. This theorem de- 
pends strongly on the spatial homogeneity of the Bianchi space-times and the 
fact that the perfect fluid flow lines can be written as conformal geodesies. 
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It remains to be seen how one could derive such theorems for other matter 
models or for inhomogeneous space-times. 
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